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Low-temperature measurements of the ac magnetic susceptibility along the [110] direction of a 
single crystal of the dipolar spin ice material Ho2Ti207, in zero static field, are presented. While 
behavior that is qualitatively consistent with previous work on Ho2Ti2 07 and the related material 
Dy2Ti2 07 is observed, this work extends measurements to appreciably lower temperatures and 
frequencies. In the freezing regime, below 1 K, the dynamics are found to be temperature activated, 
thus well described by an Arrhenius law with an activation energy close to 6J c s, a result that is 
not easily explained with the current model of magnetic monopole excitations in dipolar spin ice. 
The form and temperature dependence of the ac susceptibility spectra are found to be nontrivial 
and distinct from standard glassy relaxation. Particular attention has been paid to correcting for 
the demagnetization effect, which is quite significant in these materials and has important, even 
qualitative, effects on the susceptibility spectra. 
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I. INTRODUCTION 

Recent years have seen a flurry of activity stemming 
from the recognition 1 - that the unusual ground-state 
correlations of dipolar spin ice should support excita- 
tions that behave as magnetic charges. These magnetic 
charges, or deconfincd monopolcs, are formed through 
the fractionalization of a single spin-flip (dipolc) exci- 
tation. By virtue of the underlying correlations of the 
dipolar spin ice 'vacuum', the monopoles experience a 
Coulomb force between them and are found to flow un- 
der the influence of the magnetic fields Several experi- 
mental works and theoretical analyses of previous exper- 
imental work 2 - - — have since verified the existence of these 
monopole excitations. Furthermore, lithographically pat- 
terned, artificial spin ice systems have been created 8 - and 
also shown to exhibit monopole defects.— ~— In particular, 
the slow relaxation that is observed in the ac magnetic 
susceptibility of spin iceiSr— has recently been attributed 
to the freezing out of dynamics of deconfined magnetic 
monopoles £ 

Here, we present detailed low temperature ac suscep- 
tibility measurements, x(/), on single crystal Ho2Ti2C>7. 
These measurements represent an important step in the 
study of monopole physics of dipolar spin ice for several 
reasons, (i) Ho2Ti20 7 has received less attention, at least 
from thermodynamic measurements, than Dy2Ti20 7 . 
Ho2Ti 2 7 ac susceptibility measurements thus far have 
been limited to temperature scans of x m higher temper- 
ature and frequency regimesi 16 ' 18 ! 19 (ii) Our experiments 
represent the first low temperature ac susceptibility mea- 
surements performed on single crystal spin ice in zero 
static field along the [110] direction, (iii) In previous 
work-^ ~ 15 i 17 on Dy2Ti2C>7, there has not been a clear ac- 
counting of the demagnetization correction, which turns 



out to be rather important for quantitatively determin- 
ing the intrinsic time constant of the system as a function 
of temperature and matching to theoretical predictions^ 
The demagnetization correction also has qualitative ef- 
fects on the shape of the frequency spectra, (iv) This 
work extends measurements to much lower frequencies 
than previous work on any variety of spin ice, therefore 
delving deep into the Coulomb phase. It is our hope that 
it will provide an important benchmark by which to test 
theories of monopole physics in dipolar spin ice. 

In the pyrochlore materials Ho2Ti20 7 (HTO) and 
Dy 2 Ti 2 7 (DTO), the magnetic ions Ho 3+ and Dy 3+ oc- 
cupy a lattice of corner sharing tetrahedra. The crys- 
tal field acting on those ions creates a strong Ising 
anisotropy, with the easy axis, known as the local [111] 
axis, pointing directly in or out of the tetrahedra. For 
ferromagnetic nearest neighbor interactions, the spins are 
highly frustrated and possess a macroscopically degener- 
ate set of ground states, consisting of two spins point- 
ing into and two spins pointing out of each tetrahe- 
dron. 20 This situation is directly analogous to the pro- 
ton bonds in water ice, hence the name 'spin ice'. In 
both spin and water ice, the macroscopic degeneracy of 
ground states gives rise to an extensive residual entropy 
So = i?/21n(3/2) ) 2i^ known as Pauling's entropy^ 
In HTO and DTO, the nearest neighbor (n.n.) ferro- 
magnetic interaction is a result of dipolc-dipolc interac- 
tions, thus these systems are often referred to as dipolar 
spin ice^ 2 . Considering the antifcrromagnctic n.n. ex- 
change interaction J and the Ising anisotropy of the 
spins, the relevant effective n.n. interaction becomes 
J c ff Sf S * where 

J cB = m 2 (5D-J)/3, (1) 
m is the magnetic moment of the Dy 3+ or Ho 3+ ions and 
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D is the strength of the n.n. dipolar interaction. The fur- 
ther neighbor dipolar interactions are very well screened 
in this system.— ^ Nonetheless, they are found, theoreti- 
cally, to select a unique ground state for the system and 
should, in principle, result in a sharp, first-order phase 
transition.— However, no such transition has been found 
in experiment,— ~— possibly because of the difficulties of 
reaching equilibrium at temperatures well below 0.5 K. 

The unusual disordered ice-like state of dipolar spin 
ice, with its two-in, two-out tetrahedra, can be thought 
of as a divergence-free V • H = vacuum, with fluctu- 
ating magnetic field lines.— The simplest excitation out 
of this vacuum is the flipping of one spin. This single 
spin flip affects two neighboring tetrahedra resulting in 
one with 3-in, 1-out and the other with 1-in, 3-out. This 
removes the divergence-free character of the system and 
leads to two magnetic charges, or monopoles, of opposite 
sign, centered on those tetrahedra. In the simple nearest 
neighbor model of spin ice, these monopole excitations 
are, once created, free to roam about on the lattice with- 
out any energy cost. In dipolar spin ice, however, the 
longer range interactions give rise to a Coulomb interac- 
tion, ~ 1/r, between magnetic monopolcsi&i 

This has been associated recently^ with the slow re- 
laxation observed in the ac magnetic susceptibility of spin 
iceJ^r— At temperatures ranging from about 2 K to 6 
K, where the monopole density is high but there are very 
few double defects (4-out or 4- in) , DTO ac susceptibility 
data from Ref. [3 is well parametrized by an Arrhenius 
law with 

r = T cxp{E A /T) = r cxp(2 J ofI /r), (2) 

where 2 J ff is the energy cost of a single monopole defect. 
As the temperature is lowered and the concentration of 
monopoles decreases, the relaxation becomes slower than 
predicted by Eq.[2] A reasonably good fit of experimental 
data is obtained by performing Monte Carlo simulations 
of a gas of monopole charges on the diamond lattice in 
the grand canonical ensemble with a largely temperature 
independent chemical potential 

Despite this strong verification of the monopole hy- 
pothesis, there remains some discrepancy between the- 
ory and experiment at lower temperatures (below ~ 1 
K) in that the experiments find slower relaxation than 
is predicted by theoryji It is, thus, important to further 
explore the low temperature freezing of spin ice systems. 
Ensuring an accurate correction of the demagnetization 
effect and taking measurements to lower frequencies, as 
performed here, will permit careful comparison of exper- 
imental results with the theory of monopole defects. In 
contrast to the predictions of theory, our results show a 
temperature activated regime where the dynamics are de- 
scribed well by a surprisingly simple Arrhenius law with 
6 Jeff activation energy at low temperatures. This study, 
performed on Ho2Ti2 0y, as opposed to the more com- 
monly studied Dy2Ti207, also permits the observation 
of monopole physics with a different set of Hamiltonian 
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FIG. 1: (color online) Frequency scans of x' an d x" 01 
Ho2Ti207 with the magnetic field aligned along the [110] crys- 
tal orientation. The data shown here has been corrected for 
demagnetization and was taken on a needle shaped sample 
(sample 1). Slow relaxation is observed, characterized by a 
suppression of x a t higher frequencies and a broad peak in 
X" 



parameters, namely the spin flip rate, effective n.n. in- 
teraction, J c ff, and the monopole charge, Q. 

II. EXPERIMENT 

The single crystal Ho2Ti207 samples studied here were 
prepared at McMastcr University. The samples are from 
the same crystal growth emplo yed in the neutron scat- 
tering studies reported in Ref. l32t They are grown by 
a floating-zone image furnace technique, the details of 
which are provided in Ref. l33l . Two sample geometries 
were measured, with dimensions 1.1 x 1.1 x 2.6 mm 3 and 
0.6 x 0.6 x 3 mm 3 respectively. In each case, the longest 
side of the crystal was the [110] crystal direction, also the 
direction of the applied ac magnetic field. 

Susceptibility, x, measurements were performed using 
a magnetometer based on a superconducting quantum 
interference device (SQUID), mounted on a dilution re- 
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frigerator. The second-order gradiometer sensing coils 
are contained in the bore of an excitation coil consist- 
ing of 375 turns of NbTi wire wound on a phenolic form 
with which an applied ac magnetic field, not exceeding 
20 mOe, is generated. The gradiometer is coupled to 
the SQUID by means of a superconducting flux trans- 
former. Further refinement of the magnetometer bal- 
ance is achieved with a trim coil coupled to one branch 
of the gradiometer, in parallel with the excitation coil. 
The SQUID is contained within a small lead shield, and 
the entire gradiometer within another, larger lead shield. 
Further shielding is provided by a cryogenic /x-metal 
shield surrounding the vacuum can of the cryostat and 
by a room-temperature /i-metal shield surrounding the 
liquid helium dewar. The SQUID and SQUID controller, 
with 100 kHz modulation frequency, were obtained from 
the company ezSQ UIDM- A lock-in amplifier provides an 
ac source, and reads the feedback output from the SQUID 
controller as the resulting signal. 

The samples were mounted on a sapphire rod using 
General Electric varnish, the long edge aligned with the 
direction of the applied magnetic field with an estimated 
accuracy of ±2 degrees. The sapphire rod is clamped 
into the copper base of the sample holder and this sample 
holder is in turn heat sunk to the mixing chamber of a 
dilution refrigerator. 

Frequency scans were taken by controlling the cryostat 
at a given temperature in the range of 500 to 1300 mK, 
and sweeping the frequency in the range of 1 mHz to 
500 Hz. It was ensured that, for all measurements pre- 
sented here, the apparatus and samples were in thermal 
equilibrium. Reproducibility, and the absence of further 
thermal relaxation, was verified by taking multiple scans 
separated by several hours. The lower temperature limit 
of our data corresponds roughly to the point at which 
the maximum in \" reaches our lowest measurement fre- 
quency of 1 mHz. This temperature limit was chosen 
since much less information will be gained below that 
point without employing prohibitively low frequencies 
of measurement. No problems of thermal equilibration 
were encountered that would have otherwise increased 
the base temperature of our measurement. The exci- 
tation power was also varied by more than an order of 
magnitude to rule out heating of the sample or appara- 
tus or other nonlinear effects. Temperature scans of \ 
at fixed frequency were also obtained for four different 
measurement frequencies: 0.1, 1.2, 10 and 40 Hz. 
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FIG. 2: (color online) Superimposed absorption spectra ob- 
tained by plotting x"/XMax against ///Max- (a) The absorp- 
tion spectra taken on the less elongated sample (sample 2) 
of dimensions 1.1 x 1.1 x 2.6 mm 3 , before correcting for the 
demagnetization effect. The raw data, plotted in this way, ex- 
hibits a narrowing of the absorption spectrum with decreasing 
temperature, (b) Data fully corrected for the demagnetization 
effect, taken on a needle-shaped sample (sample 1). The spec- 
tra change qualitatively, with a clear broadening of the low- 
frequency tail as the temperature is reduced. However, the 
high-frequency tail continues to narrow with reducing tem- 
perature. Above 1 K, the low-frequency tail also shows slight 
narrowing with reducing temperature, although the effect is 
rather subtle. Arrows indicate increasing temperature. 



III. RESULTS 

Results of ac susceptibility frequency scans, shown in 
Fig. [1] are qualitatively consistent with previous results 
taken on the spin ice material Dy2Ti20yF^r— Slow re- 
laxation is observed, characterized by broad absorption 
spectra x"if) an d suppression or blocking of the in phase 
susceptibility x' at higher frequencies. Rapid freezing of 
the magnetic moments can be seen through the sharply 



dropping peak frequency, /Max, of the absorption spec- 
trum x"(f) as the temperature is lowered. We see that 
by 500 mK, /Max is already below our frequency window 
at less than 1 mHz. At higher temperatures approaching 
1.3 K, the peak position is beginning to plateau at around 
100 Hz. The magnitude of the susceptibility monotoni- 
cally increases with decreasing temperature. However, 
judging by the peak height of x"> the susceptibility ap- 
pears to be leveling off at the lowest temperatures stud- 
ied. 
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FIG. 3: (color online) Frequency versus inverse temperature 
taken from this work and fits with different methods. Green 
diamonds are obtained from the maxima Tmelx(/) of temper- 
ature scans of x" ■ The corresponding best fit Arrhenius law 
is shown as the green line. The blue circles are obtained from 
the maxima of the absorption spectra, /Max(T). The best fit 
Arrhenius law, shown in magenta, gives EaI^b = 10.70 K, 
which is very close to 6J e fi for Ho2Ti207. In fact, an Arrhe- 
nius law with EaI^b = 6J e a also fits the data to within the 
uncertainty, as shown in blue. The fits are plotted as solid 
lines in the range of the fitted data and extended with dashed 
lines. The 2J c g Arrhenius law that should be expected in the 
plateau region is shown in red. 



In order to make comparison of the spectra more 
amenable, we superimpose the spectra in Fig. [2l This 
is done by plotting data normalized by the peak sus- 
ceptibility and frequency, so plotting x"/XMax against 
///Max- Interestingly, if this is done before correcting 
for demagnetization with the less elongated sample, one 
finds a subtle narrowing of the absorption spectrum with 
lower temperatures, as shown in Fig. [2ja). This effect 
was previously noticed by Snyder et aL— in measuring 
the susceptibility of Dy 2 Ti 2 0y. 

Because of a large magnetic moment on the Ho 3+ sites, 
the susceptibility of Ho 2 Ti 2 07 is quite large and the de- 
magnetization effect is crucially important to consider. 
In fact, without correcting for demagnetization, one can 
obtain qualitatively different spectra due to mixing of x' 
and x" ■ If a demagnetization correction is performed, 
the narrowing behavior becomes less obvious, as seen in 
Fig. [2{b). The low-frequency tails of x"(f) are found to 
broaden with decreasing temperature, whereas the high- 
frequency tails still narrow significantly with decreasing 
T. Even in the demagnetization corrected data, there 
is some hint of the low-frequency side of the spectra be- 
ginning to broaden very slightly above 1 K. Thus, the 
relaxation in spin ice is found to be noticeably different 
from the behavior of a standard glass or spin glass, where 



the absorption spectra show clear and largely symmetric 
broadening with reduced temperaturei 35 ' 36 

Quantitatively, the broadening of the absorption spec- 
tra on a log scale can be parametrized by the half 
width at half maximum, either on the high frequency 
side (HWHM+) or low frequency side (HWHM_), both 
of which are plotted against temperature in the inset 
of Figure [3] HWHIvL is roughly constant at around 
0.9 decades from 1.3 K down to 0.85 K (though show- 
ing barely discernible narrowing with lower T). Below 
that point, it shows broadening, to almost 1.2 decades 
by 650 mK. HWHM + shows a stronger and opposite 
temperature dependence, increasing from ~ 1 decade at 
the lowest temperatures studied here to 1.5 decades at 
1.0 K. These numbers can be compared to the 0.7 decades 
HWHM and symmetric absorption spectrum that results 
from a single energy barrier to relaxation. The asym- 
metry of the spectra indicate that the relaxation is not 
described by a single characteristic Debye form. Simple 
empirical forms such as the Cole- Cole j2I Davidson- Cole 3 ^ 
and Havriliak-Negami 3 ^ functions are also not good fits 
to the data at low T, unlike what was found above 1.8 K 
in a previous analysis of the susceptibility of Dy 2 Ti 2 0y J£ 

At low temperatures, the peak absorption frequency 
appears to approach an Arrhenius law, /Max = 
/o exp(— Ea/T), thus resulting in a straight line when 
log /Max is plotted against l/T, as seen in Fig. [3] Fitting 
the data below 1 K gives a best fit Arrhenius law with 
activation energy E A = 10.70 K and f Q = 5.74 x 10 5 
Hz. The resulting reduced chi-square statistic, X 2 , is 
0.12, implying that it is a very good fit to the data. 
The value of Ea so obtained is tantalizingly close to 
6J off = 10.80 K for Ho 2 Ti 2 7 3 In fact, an Arrhenius 
law with exactly Ea = 6J c s is also a more than adequate 
fit with f Q = 6.57 x 10 5 Hz and reduced X 2 = 0.18. The 
low values of X 2 suggest that the fit is somewhat under- 
constrained. Taking a 95% confidence interval allows for 
Ea = 10.70 ± 0.15 K. Above 1 K, the temperature de- 
pendence of the frequency /Max becomes shallower than 
the Arrhenius-law fit, likely moving towards the quasi- 
plateau regime that was seen in DTO . 12 ' 14 At the lowest 
measured temperatures, an onset of Arrhenius behavior 
was seen in DTOi^ and in a later work a ~ 6J e g Ar- 
rhenius law was shown to fit the lower end of that same 
experimental data.— However, insufficient low tempera- 
ture data has been obtained on DTO to be reasonably 
certain that the Arrhenius behavior develops as is seen 
here in HTO. 

Detailed temperature scans at four different frequen- 
cies of measurement, 0.1 Hz, 1.2 Hz, 10 Hz, and 48 Hz, 
are shown in Fig. |4j Once again, the hallmarks of slow re- 
laxation are observed, with suppression of x' at l° w tem- 
peratures and a peak in x"{T) near the inflection point 
of x'{T). The peak in x" shifts to higher temperatures 
and broadens as the frequency of measurement, /, is in- 
creased. The peak position of the temperature scans, or 
lMax(/), may also be extracted in order to characterize 
the relevant time scale as a function of temperature. 
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FIG. 4: (color online) Temperature scans of x' an d x" °f 
Ho2Ti2 07 with the probe field aligned along [110], for four 
different measurement frequencies: 0.1 Hz, 1.2 Hz, 10 Hz and 
48 Hz. This data has been taken on a needle-shaped sample 
(sample 1) and has been corrected for the demagnetization 
effect. Again, slow relaxation is observed, with a maxima in 
x' or x" moving to lower temperatures with lower frequency 
of measurement. 



It is found that the peak positions of the frequency 
scans and of the temperature scans lead to different func- 
tional forms, so that fMax(T) ^ /(Tkiax)- This is most 
easily demonstrated with surface plots of T) and 
X"(/j T), shown in Fig. [5] The maxima in %" taken along 
the frequency axis and along the temperature axis are 
shown in unfilled and filled circles respectively and do 
not sit on the same curve. While there is a very signifi- 
cant difference between these paramctrizations at higher 
temperature, the two curves appear to be tending to- 
ward a single function at lower temperatures. In previous 
workj£ /Max(r) was used to compare with theory. Ide- 
ally one would like to have theoretical calculations that 
can reproduce the entire %(/, T) surfaces that are shown 
in Figure [H but for now, fMax(T) will serve as a good 
representative parameterization of our data. 

The demagnetization correction shifts the peak fre- 
quency significantly lower, especially at lower tempera- 
tures where the static susceptibility becomes quite large. 
In other words, the demagnetization effect, increases the 
measured relaxation times relative to the relaxation times 
that would be measured from the true, bulk susceptibil- 
ity. This is a result of x' feeding into the measurement 
of x" ■ This reduces the apparent absorption spectrum 



x'a a t low frequencies where x' is iar g e ; thereby pushing 
the peak frequency /Max to higher frequency. For more 
details on the demagnetization effect and correction, see 
the Appendix. 



IV. COMPARISON WITH PREVIOUS WORK 

Multiple experimental techniques have been used 
to study the relaxation of the three well established 
spin ice materials: Dy 2 Ti 2 7 (DTO)^^^ Ho 2 Ti 2 7 
(HTO)i&±&2^ and Ho 2 Sn 2 7 (HSO)J£ In ac suscep- 
tibility, neutron spin echo, neutron scattering, and /iSR 
relaxation, the following general features are observed. 

Above ~ 15 K, one observes a high temperature regime 
where the transition rate of the individual spins is chang- 
ing rapidly in T. The energy scale for this regime ap- 
pears to be set by the energy gap to the next excited 
crystal field energy level at ~ 210 K for DTQi^r— and 
~ 290 K for HT041 The Ho 3+ (or Dy 3+ ) moments have 
a truly Ising doublet ground state with matrix elements 
of the local J x and J y angular moment operators iden- 
tically zero. Transitions within this doublet are thus 
forbidden transitions. Spin flips must therefore occur 
through either quantum tunneling (a slow process) or 
by passing through the excited crystal field energy levels 
(a much faster process). Relaxation in the high tem- 
perature regime can thus be described by an Arrhenius 
law with an activation energy of ~ 300 K. As the tem- 
perature is lowered it becomes difficult to populate the 
excited crystal field states reducing the transition rate. 

Below ~ 15 K, the spins are very Ising-like and the 
materials enter a quasi-platcau regime where the tun- 
neling rate levels off. The vast majority of spin flips 
are now occurring via tunneling processes. Note how- 
ever, the tunneling is also accompanied via a change in 
energy resulting from monopole excitations, thus there 
remains a temperature dependence to the relaxation. 
This plateau should correspond to the high temperature 
tail of a 2 J G ff Arrhenius behavior that results from the 
thermal excitation of independent (deconfined) monopole 
defectsi&i The monopole density in the plateau regime 
is high enough that there is significant screening of the 
magnetic charges to make the nearest neighbor spin ice 
model largely applicable. 

Approaching J c g (around 1.5 to 2 K)^ the system en- 
ters a regime exhibiting rapid freezing out of dynamics^ 
Here the time scales begin to increase dramatically as the 
system condenses into a disordered 2-in-2-out configura- 
tion. Equivalently, it is now becoming thermally difficult 
to form monopole defects. The monopole density de- 
creases which in turn reduces screening of the magnetic 
charges, making full treatment of the dipolar spin ice 
model necessary to describe the dynamicsi&£ 

While different techniques and the different spin ice 
materials give quantitatively different time scales, they 
all agree on the qualitative behavior and on three clear 
regimes of relaxation. For example, HTO neutron scat- 
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FIG. 5: (color online) Surface plot of \ (top) and x" (bottom) 
with temperature T on the y-axis and measurement frequency 
/ on the a;-axis. Also shown, are the positions of the max- 
ima in x" either taken from temperature scans (black, filled 
circles) or frequency scans (red, open circles). 



tcring experiments probing particular wave vectors re- 
solve fluctuation rates ~ 10 ns, that are several orders 
of magnitude faster than those seen in the bulk suscep- 
tibility, with time scales around ~ 1 msj22 Despite this 
significant difference in relaxation times, the same quali- 
tative shape of f(T) described above is observed^ Neu- 
tron scattering measurements by Zhou et al. also found a 
quantum tunnelling regime on the time scale of ~ 0.01 ns 
in the dynamic spin ice material Pr2Sn207i^ Ehlers et 
al. have found, using neutron spin echo experiments, 
the same three regimes of relaxation ! 18 ' 41 Susceptibility 
measurements of DTO— have found the time scale of 
the plateau to be on the order of ~ 1 ms, while /xSR 
measurements of DTO have found the above qualitative 
shape, but with the time scale of the plateau relaxation 
on the order of 0.1 [is£H This mismatch of timcscalcs 
in the dynamics as determined by different measurement 
techniques remains an open question in the study of spin 
icei 28 ' 40 The answer may lie in the wavevectors that are 
accessed by these various measurements; while neutron 
scattering accesses specific wavevectors, bulk susceptibil- 
ity addresses k = and /iSR is a local measurement, thus 
has contributions from all wavevectors^ 

Interestingly, measurements of DTO by Orendac et al. 
using the magnetocaloric effect have accessed a frequency 
and temperature range lower than that of any ac suscepti- 
bility measurements^ Those measurements also seem to 
access a different time scale, showing both a higher rate 



of relaxation and a shallower temperature dependence of 
the freezing than what is observed in ac susceptibility. 
Henceforth, we concentrate on comparing our data with 
that of other bulk susceptibility measurements of spin 
ice. 

Even between ac susceptibility measurements, a direct 
comparison of data can be difficult to achieve. For many 
of the measurements published on spin ice, a correction 
of the demagnetization effect has not been performed, 
or is not incorporated into the presented data. Many of 
the previous results were obtained on powder samples, 
with which it is difficult to compare aligned single crys- 
tal measurements. Of course different spin ice materials 
have different parameters in their respective Hamiltoni- 
ans, which makes their quantitative behavior different. 

Finally, comparison is further hindered by the com- 
plex dynamical behavior of spin ice and the various 
measurement schemes and parametrizations that can be 
used. One of the most obvious ways to parametrize 
the data is with the maximum of the absorption spec- 
tra at fixed temperature, or fMax(T). This parameter 
has been employed in several experimental works and in 
theory i 6 ' 7 ' 14 ' 15 Alternatively, temperature scans at con- 
stant frequency may be employed and two maxima in x" 
can be resolved at T Low (/) and T mgh (f)~£*£ The higher 
temperature maximum, Tnigh, appears to be present only 
for frequencies above a certain threshold (near 200 Hz 
for DTQi 5 -). While T Hig h has not been observed in HTO 
in zero field due to an unsuitable range of temperature 
and frequency, it has been observed, by Ehlers et a/] 18 i 41 
under an applied field of 1 T. In temperatures scans, 
one may also extract maxima in x' I n Figure [HI we 
show a compilation of several of the previous results on 
DTO^ii^i 5 . HTOi 6 - and HSOi 6 . compared with our re- 
sults on HTO. 

While qualitatively similar behavior is observed in ac 
susceptibility for all the materials studied, quantitative 
differences remain. There is very little difference in the 
data below 1 K between our results on HTO and re- 
sults on DTO from Refs. [T^ITfjl . despite the differences 
in Hamiltonian, specifically the effective n. n. exchange 
energy scale J e g, and in demagnetization effects. We can 
draw from the physics of monopole excitations discussed 
in Refs. QUO] to suggest what relations between the dif- 
ferent relaxation rates should be expected. 

First, the tunneling rates will be different in the two 
materials studied. Tunneling rates for Ising moments like 
those in spin ice arc functions of several different param- 
eters. The symmetry and strength of the crystal field, 
hence the energy of excited states above the Ising dou- 
blet, play an important role, giving the energy barrier 
through which the moments must tunneL 12 ' 14 Compo- 
nents of the exchange and dipolar interactions transverse 
to the local Ising direction result in mixing with the 
excited crystal field energy levels, permitting tunneling 
to occuri 36 i 45 Furthermore, nuclear hypcrfine interaction 
strengths, very different between Dy 3+ and Ho 3+ , can 
play a role in assisting or blocking spin ffipsi 36 ' 45 To fur- 
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ther complicate the system, we may anticipate that there 
will not be a single microscopic tunneling time, but a dis- 
tribution. However, spin ice has so far been treated with 
a single spin flip rate fo&2- In the regime where tunnel- 
ing is dominant (from temperatures in the quasi-plateau 
regime and below), increasing the spin flip rate would 
be equivalent to stretching the frequency axis. It is not 
clear at this time what the difference in tunneling rates 
is between the various spin ice materials, though our re- 
sults suggest that they are at least in the same order of 
magnitude. 

Within the nearest-neighbor spin ice (NNSI) model, 
the only other parameter governing the dynamics of the 
system would be J c s, or the energy barrier that must be 
overcome in order to excite monopole defects.— Chang- 
ing J e fj would stretch the temperature axis, either elon- 
gating or shortening the plateau regime. In Dy2Ti207, 
D nn = 2.35 K and J nn = -1.24 K giving J cS = 1.11 KM 
In Ho2Ti2 07, the dipolar interaction is essentially the 
same, whereas J nn = —0.52 K, giving a much higher 
value of J ff = 1.83 Because of the higher J e g in 
HTO, it enters the freezing, 2-in-2-out regime, at higher 
temperatures than DTO<^ 

When one goes beyond the NNSI model into the dipo- 
lar spin ice model, the monopole defects begin to inter- 
act with a Coulomb law, V(ry-) = (fJft/ArfQiQj/rij, as 
long as they occupy distant sites on the diamond lattice. 
Thus a third parameter, the magnetic charge Q, starts 
to play an important role at lower temperatures as the 
monopole density is reduced and there is less screening of 
the magnetic charges. The monopole charge is not that 
of the conjectured monopolcs of particle physics, but is 
much smaller and has been determined theoretically^ and 
experimentally 2 - to be 
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V 3/^0 
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where fx = 2cjji±b{Jz) is the magnetic moment of the 
rare-earth ions, ad is the distance between adjacent sites 
of the diamond lattice and a is the distance between spins 
on the pyrochlore lattice. Q is not directly related to J e g, 
but only the strength of the dipolar interaction D nn . If 
the NNSI model were to hold, a plot of /Max To versus 
T / J e f{ would in principle give the same result for all spin 
ice samples, from the quasi-plateau regime and below. In- 
troduction of this third parameter, the monopole charge 
Q, should make such a simple scaling of relaxation curves 
impossible. 

Surprisingly, despite all the differences between sam- 
ples and measurements, our frequency scan maxima /Max 
are quite close to those of Snyder et a/. 12 ' 14 on a poycrys- 
tallinc sample of DTO. While overlapping temperatures 
are not available, our results do seem to be approaching 
those of Matsuhira et al^, also on polycrystalline DTO, 
at higher temperatures. 

We may more directly compare our temperature scans 
with the results of Ref . [la where HTO has also been stud- 
ied and a demagnetization correction has been applied. 



Shown in the inset of Figure [6l are our results at / = 10 
Hz on a single crystal oriented along [110] and those of 
Matsuhira et aL— at the same frequency on a powder 
sample of HTO. There is a clear difference between the 
curves, suggesting that the magnetic field orientation is 
an important parameter. The relaxation times deter- 
mined from Tlow(/) are found to have a steeper tem- 
perature dependence for the polycrystalline sample than 
for the [110]-oriented sample. This would seem to sug- 
gest that energy barriers to the movement of monopole 
defects may be lower along the [110] and symmetry re- 
lated directions (in other words, along the so-called^ a 
or P chains) than in other crystal directions. In a perfect 
gas of interacting charges all directions are equivalent, 
but the magnetic monopoles in spin ice are confined to 
travel only in certain paths&£ which could give rise to 
anisotropic behavior. 

Shi et aZ.jii have studied the relaxation of two DTO 
single crystals, one oriented along the [111] direction and 
the other oriented along the [100] direction. They note 
that there is a difference between Tf(/) (the maximum 
in x'(T)) between these two scenarios. However, it is not 
clear that their data has been corrected for demagneti- 
zation. The two samples have different geometries and 
the different orientations will give rise to different mag- 
nitudes of susceptibility since the magnetic field will be 
differently aligned with the different basis spins. This 
change in magnitude will heavily influence the apparent 
susceptibility xaU, T). 



V. CONCLUSIONS 

To conclude, we have performed a careful study of 
the frequency dependent spin freezing in a single crys- 
tal of the dipolar spin ice material Ho2Ti2 07 at low 
temperatures. The results are qualitatively consistent 
with research by other groups on Ho2Ti207,— on the 
related materials T)y 2 Ti 2 7 ^^. and Ho 2 Sn 2 07^ and 
also with theoretical work that relates such dynamics 
to the motion of magnetic monopole excitations The 
results presented here represent an exploration much 
deeper into the Coulomb phase of spin ice, with much 
lower temperatures and frequencies of measurement than 
were employed for ac susceptibility measurements previ- 
ously. 

One of the most striking aspects of our results is the 
observation of a seemingly temperature activated regime 
below 1 K where the dynamics are well fit by an Arrhe- 
nius law, consistent with a 6 J e ff activation energy. It has 
also been noted^ that an Arrhenius law with an activa- 
tion energy close to 6 J e ff is consistent with the low tem- 
perature limit of measurements on Dy2Ti207 . 12,14 While 
very few points contribute to that conclusion for DTO, it 
nonetheless supports our observation of such an effect in 
HTO. A more thorough investigation of the low temper- 
ature behavior of DTO with full demagnetization correc- 
tion would be valuable to verify or disprove the existence 
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FIG. 6: (color online) Comparison of various ac susceptibility 
data on spin ice compounds. Three methods for determin- 
ing f(T) are used. Two different maxima may be determined 
from temperature scans of x"> which are labelled as TLow and 
Thigh- Only TLow has been determined in our data because 
of the higher frequencies and temperatures required to mea- 
sure Thigh- Data on Dy2Ti2C>7 from Snyder et al^ and from 
Matsuhira et a/.— are shown, and discrepancies can likely be 
attributed to demagnetization effects. Data from Matsuhira 
et a/.— are shown on the materials H02T12O7 and Ho2Sn2 07 
as well as our results on Ho2Ti2C>7. The inset shows temper- 
ature scans (scaled arbitrarily for ease of comparison) from 
this work (red triangles) and from Ref. flrl at the same fre- 
quency of 10 Hz. The discrepancy could be attributed to a 
difference between powder samples^ and the aligned single 
crystals measured here. 



tion for such materials and its impact on the shape of 
the absorption spectra, for example. 

When compared to experiments on Dj^T^Oy, theory 
currently shows a persistent mismatch at lower temper- 
atures, with the experiments showing longer relaxation 
times than would be expected theoretically.— While those 
experiments have not taken into account the demagneti- 
zation effect , this will not help explain the discrepancy 
and should in fact make it worse as correcting for demag- 
netization will lower the peak frequencies. 

In order to resolve such discrepancies, it has been 
suggested^ that higher order terms, such as next near- 
est neighbor exchange, J2, may need to be added to 
the currently accepted Hamiltonian of the dipolar spin 
ices. Such interactions were already found to be quite 
important for matching theory to experiment in DTO. 46 
It would also be useful to consider the possibility that 
the low temperature susceptibility of Ho2Ti207 is af- 
fected by bound pairs of monopolcs. Such very slowly 
relaxing pairs of monopoles have recently been predicted 
to occur after thermally quenching spin ice4£ It is not 
clear, however, whether experimentally inaccessible cool- 
ing rates are required to form these pairs or whether they 
are practically unavoidable in the real systems. 

It is hoped that the very exciting theory of monopole 
excitations in spin ice can begin to explain some of the 
peculiar effects that are seen in these susceptibility mea- 
surements. Equivalcntly, these results, having been per- 
formed with careful demagnetization correction on single 
crystal samples and to comparatively low frequencies and 
temperatures, should serve as an important benchmark 
for testing the theory of monopole physics in spin ice. 



of the simple law that we have observed. 

It is not clear what could give rise to such a simple 
temperature activated regime in dipolar spin ice, which 
seems at odds with the current theoretical treatment 
of magnetic monopoles in dipolar spin ice.— i Such the- 
ory instead predicts ever increasing activation energies 
as monopole density and thus screening is reducedi&i 
Furthermore, a one-parameter ( J e ff) fit to this low tem- 
perature regime surprisingly suggests that the monopole 
charge Q might be an irrelevant parameter. There re- 
mains the possibility that the accessible temperature 
range is simply narrow enough that a more complex 
/(1/T) function only appears to follow a linear behav- 
ior. The Arrhcnius law holds well over almost 4 orders 
of magnitude in frequency, though applies only over a 
factor of 2 in temperature. 

Our measurements have also illustrated the complex 
nature of the spectra. While qualitatively glassy relax- 
ation is seen, the precise behavior of the susceptibility 
spectra is quite different from that of spin glasseSi 35 ' 36 
In particular, as the temperature is lowered, the high 
frequency tail is found to narrow appreciably where the 
low frequency tail broadens. Clearly illustrated here is 
the extreme importance of the demagnetization correc- 
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Appendix: Calibration and Demagnetization 
Correction 

Because of the large susceptibility of Ho2Ti207, the de- 
magnetization correction is quite significant and must be 
dealt with quite carefully in order to obtain accurate re- 
sults. Not accounting for the demagnetization effect can 
even result in qualitatively different susceptibility spec- 
tra. Here we use the standard relation 

- = — -AnN (AT) 
X XA 
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FIG. 7: (color online) Comparison of the calibrated and de- 
magnetization corrected susceptibility of the two sample ge- 
ometries studied in this work. The needle-shaped sample with 
demagnetization factor JVi is shown with a dashed red line 
and the less elongated sample with demagnetization factor 
N2 is shown as a solid blue line. Excellent agreement be- 
tween the two data sets is obtained after suitable selection of 
the calibration factors. While overlap is determined from the 
absorption spectra x" > the in-phase susceptibility \' is also 
found to match very well. 



to perform the demagnetization correction, where \ is 
the true complex susceptibility of the sample, \A is the 
apparent, or measured susceptibility, and N is the de- 
magnetization factor. Because x = x' ~ i s complex 
and has a great deal of frequency dependence, there is 
significant mixing between x! an d x" ■ This alters the 
shape of the spectra obtained. Thus it can be seen from 
expanding Eq. IA.11 that 



X 



X 



x' A -^N(x2 + xf) 



(1 



x'k 



(1 - AttNx'a) 2 



(AttNx'X) 2 



(A.2) 



(A.3) 



as discussed in Ref. [l5|, for example. 

While one method of calibration could involve the use 
of a calibration standard, such as a perfectly diamagnctic 
superconducting sample, this method is fraught with dif- 
ficulties. To be completely accurate, one would need a 
calibration standard that is exactly the same size and 
shape and positioned in the same place as the sample. 

Our method, instead, involves the measurement of two 
different sample geometries with demagnetization factors 
Ni and N 2 respectively. The raw susceptibility signal 
of each sample geometry is taken in arbitrary units, as 
V\ and V2 respectively, over a common range of tem- 
peratures and frequencies. An array of trial calibration 
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FIG. 8: (color online) The effects of the demagnetization ef- 
fect on the susceptibility of Ho2Ti2 07 at two different temper- 
atures, 600 and 800 mK. Blue and red curves are uncorrected 
data (apparent susceptibility) measured on a short sample 
of dimensions 1.1 x 1.1 x 2.6 mm 3 and on a needle shaped 
sample of dimensions 0.6 x 0.6 x 3 mm 3 , respectively. The 
black curves are data that has been fully corrected for the 
demagnetization effect. As can be seen from the figure, the 
correction is very significant, more than doubling the suscep- 
tibility at low temperatures and low frequencies where it is 
most important. The inset shows the maximum frequencies 
/Max obtained before (blue and red) and after (black) correct- 
ing for demagnetization, as a function of inverse temperature. 



factors C\ and C2 for each of the samples is generated. 
For each entry, the apparent susceptibilities, xai = G{V\ 
and XA2 = C2V2, are corrected for demagnetization ac- 
cording to Eq. IA.ll An overlap is determined from 
the absorption spectrum data points with the parame- 
ter XXx'i ~ x'i) 2 / (x'{ 2 + X2 2 ) an d the calibration factors 
C\ and C2 are chosen in order to minimize that parame- 
ter (maximize the overlap) or so xi = X2 = X f° r a H data 
points. 

This method utilizes a calibration of the sample which 
ensures proper scaling of the apparent susceptibility with 
sample geometry, thereby providing a very accurate cor- 
rection of the demagnetization effect. The only disad- 
vantage, besides increased measurement time, is that it 
entails the assumption that the true susceptibility of the 
samples \ is n °t sample geometry dependent and that 
only the apparent susceptibility xa changes. It is likely 
that, apart from extreme situations like a thin film of 
material, the shape of the sample does not affect intrin- 
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sic thermodynamic quantities such as the true magnetic 
susceptibility. 

Specifically, we have used, for our calibration, fre- 
quency scans of the susceptibility at three different tem- 
peratures: 600, 700 and 800 mK. The most needle-shaped 
sample geometry has dimensions 0.6 x 0.6 x 3 mm 3 and 
therefore has a demagnetization factor of N\ — 0.085. 
The other, less elongated but larger sample has dimen- 
sions 1.1 x 1.1 x 2.6 mm 3 with demagnetization factor 
N 2 = 0.171. The demagnetization correction was cal- 
culated with the analytical form for a rectangular prism 
given in Ref. H|| The excellent overlap so-obtained is 
shown in Figure [7] The effects of sample geometry on 
the apparent susceptibility spectra and on the peak fre- 
quencies, /Max, compared with the fully corrected data, 
are shown in Figure [5] 

The importance of the demagnetization effect has not 
always been highlighted in other measurements of spin 



ice though in some cases, it has been pointed out to be 
quite important^^ In measurements of polycrystalline 
HTO, by Matsuhira et al, a demagnetization correction 
was performed^ While results of Matsuhira et al. on 
DTO were corrected for demagnetization before the fi- 
nal analysis, much of the published data with which we 
might compare our results are left uncorrected.— Re- 
cent magnetization measurements 3 - were also corrected 
for demagnetization, particularly important since they 
were performed on a complicated sample geometry in 
large magnetic field. In some works, no mention was 
made as to whether a demagnetization correction was 
used^ ~ 14 ' 17 Here we have measured a single crystal with 
a well-defined and calculable^ demagnetization factor N 
and have performed a careful correction of the demagne- 
tization effect. We can thus say with confidence that we 
are measuring the true, bulk susceptibility of Ho2Ti207. 
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